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We study the semiclassical fluctuation problem around bounce solutions for a self-interacting 
scalar field in curved space. As in flat space, the fluctuation problem separates into partial waves 
labeled by an integer I, and we determine the large I behavior of the fluctuation determinants, a 
quantity needed to define a finite fluctuation prefactor. We also show that while the Coleman- 
De Luccia bounce solution has a single negative mode in the I = sector, the oscillating bounce 
solutions also have negative modes in partial waves higher than the s-wave, further evidence that 
they are not directly related to quantum tunneling. 
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The problem of false vacuum decay in the presence of gravity []] provides an important window into the behavior 
of interacting quantum fields in curved space-time, and is also important for our understanding of string theory and 
quantum gravity 0, 0, 3 , and inflationary models of cosmology |5[ . Since the pioneering work of Coleman and De 
Luccia Q , much has been learned regarding the existence and properties of bounce solutions for interacting scalar 
I> ' fields coupled to gravity 0, 0, 13, 13, lM EO, E3, E3, 

ITU Hil ITfij , and consequently the exponential factor in the false 
vacuum decay rate. On the other hand, relatively little is known about the prefactor in the decay rate. This is in 
distinct contrast to the flat space case where the entire computation is well understood physically and mathematically; 
analytically in the thin- wall limit [H [H E [H Hj , and numerically in general [U Hi IHHi] . Here we begin to 
address this prefactor question with coupling to gravity by studying the problem of quantum fluctuations about the 
£ — ' bounce solutions. A full solution to this problem is not possible at present for the simple reason that computing 
the renormalizcd fluctuation prefactor would require an understanding of the renormalization of quantum gravity. 
■ However, we argue that certain interesting things can be learned, in particular in the limit where the gravitational 
I background is fixed to be de Sitter space. 

In the flat space false vacuum decay problem, the fluctuation operator separates into partial waves labeled by an 
integer I, and there are three important types of modes. In the / = sector there is a single negative mode and 
this is responsible for the decaying nature of the problem 0, 0, ITsl ITflj . In the I = 1 sector there are four zero 
modes corresponding to translational invariance in four-dimensional Euclidean space, and these zero modes lead to 
Oh. collective coordinate contributions to the overall fluctuation determinant. For I > 2 the eigenvalues are all positive, 
and since for each I the fluctuation operator is a one-dimensional radial operator, one can compute the determinant 
straightforwardly using the Gel'fand-Yaglom method (described below in Section ITV|) . Formally, the determinant of 
the full fluctuation operator is a product of the determinants for all I, including degeneracies, so the large I behavior is 
crucial for defining a finite renormalizcd fluctuation determinant. In the thin- wall limit, where the energy gap between 
the true and false vacua is small, the computation can essentially be done analytically 0, ItH fl8l Il9j , and one has a 
beautiful physical picture of this process as nucleation of bubbles. Away from the thin-wall limit, the computation can 
be done by various approximate or numerical approaches |2ll \2'A |2'1 l24j . Our main motivation here is to investigate 
how the behavior of the fluctuation operator is affected by the inclusion of coupling to gravity. 

With the inclusion of gravity, Coleman and De Luccia argued that the bounce solutions are still radially 
symmetric (although this has not been rigorously proved, as it has been in flat space p^jl- Interestingly, new classes 
of bounce solutions arise, with different physical interpretations. The Coleman-De Luccia (CDL) bounce generalizes 
the flat space bounce and is presumed to be associated with quantum tunneling Q]. There also exists the Hawking- 
Moss (HM) bounce [(| which is interpreted physically in terms of a thermal transition 0, ^| . More recently it has 
been shown that there are also "oscillating bounce" solutions in which the scalar field passes over the barrier more 
than once 0, IT3L IT3 . Il5l| . As emphasized in > these oscillating bounces interpolate between the CDL and HM 
bounces, and reflect the thermal character of quantum field theory in de Sitter space. Since all these bounces are 
radial, a similar separation of the fluctuation problem into "partial waves" is possible, with the physically plausible 
assumption that such radial fluctuations dominate. But even with this separation, the fluctuation problem is still 
considerably more subtle with the inclusion of gravity, as it requires a detailed constraint analysis to disent ang le the 
physical fluctuation fields. Here we consider the scalar fluctuations in the formalism developed in [26II27LI2T 
existence of negative modes in the I = sector for these scalar fluctuations has been investigated in [2(| l27l 
We extend this fluctuation analysis in several ways by considering the behavior for higher I. We study two main 
questions: First, we investigate the large I behavior of the fluctuation determinants within each partial wave sector. 
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Wc find an explicit expression for the leading large I behavior, and a numerically accurate estimate for the sublcading 
behavior. Second, we analyze the existence of negative modes not just in the I = sector, but also for higher I, 
and show that the oscillating bounce solutions have negative modes for higher I. This is further evidence for the 
physical picture in that these oscillating bounces are not directly related to quantum tunneling, but rather reflect 
the thermal nature of quantum field theo ry i n de Sitter space. We are not able to study the I = 1 sector, as this 
fluctuation formalism does not apply here [2y, |22| , and so this requires a separate study. 

In Section ^ we review the model and the construction of bounce solutions to the classical equations of motion. 
In Section IIIII we summarize the scalar fluctuation problem to be studied. Section IIVI is devoted to the study of the 
large I behavior of the fluctuation determinant, in which we review the flat space approach. In Section we count 
the negative modes for various I for fluctuations about bounce solutions. Section IV II contains our conclusions and an 
Appendix gives the relation of our fluctuation operator to other forms considered in the literature. 



II. CLASSICAL BOUNCE SOLUTIONS 



Before discussing quantum fluctuations, we briefly review the derivation of the bounce solutions themselves. Wc 
consider the four dimensional self-interacting scalar field system with Euclidean action 



Se = J d 4 Xy/g 
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where the gravitational coupling is expressed as k = 8ir/M 2 ^. In terms of the proper time a, the metric has the form 



(2.1) 



ds 2 = da 2 + a 2 (a)dVL 
The classical Euclidean equations of motion are 
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where the overdot denotes and V'{4>) = SV /^ ■ K = 0, ±1 corresponding to flat, closed/open universes, respec- 
tively. The boundary conditions for the bounce solutions are 



0(0) = , 
a(0) = , 



<Mc max ) = o , 

a(a max ) = , 



(2.5) 



where er max is defined by the last equation: a(cr max ) = 0. In this paper we consider K = 1, which leads to the 
normalization condition 



d(0) = 1 . 

We choose the standard quartic scalar field potential 0, 0] : 
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which is sketched in Figure ^ The function / is a function of the dimensionless field ip = -: 
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The potential V(cj>) has two local minima, a false vacuum 0f v , and a true vacuum <j) tv , separated by a local maximum 
</>to P , chosen to be at </) top = 0. A crucial difference between the flat-space case and the gravitational case is that 




FIG. 1: The scalar field potential V(tp) with the parameters: f3 = 45, b = 0.25, e = 0.23. 



the overall constant Vtop in the potential is now significant, as it plays the role of a cosmological constant Q. A 
corresponding mass scale is defined as 



H ^ = \j—~\J^r ■ (2-9) 

The dimensionless parameter (3 in (|2.7|) characterizes the ratio of the barrier curvature (at </> top = 0) to -ff t 2 op 

/?^^F , (2-10) 

-"top 

and is an important quantity in determining the existence and form of bounce solutions 0, 0, 0, Il2l Il3| . Another 
useful dimensionless quantity is the ratio of the field scale v to the Planck mass M p \: 



(2-11) 



and we consider here values such that the potential is everywhere positive. 

The three critical points of V((f>) correspond to three trivial solutions to the bounce equations l|2.3|) - l|2.6[) . in which 
<t> is constant at one of these critical values 4> c such that V'(cf> c ) = 0: 

4>{o) = <f) c 

a{a) = — sm(Ha) . (2.12) 



The three solutions of this form are : (i) the false vacuum constant solution with <f> = <fif Vj and Hf v = y/ nV((j)f v )/3; 
(ii) the true vacuum constant solution with = <p tv , and H tv = \J nV((j)t v )/3; (iii) the Hawking-Moss solution 
with 4> = 0to P = 0, and H top given by (|2.9|) . 

More interesting are the bounce solutions in which <p is not constant. For definitencss, we restrict our attention to 
bounces beginning near the true vacuum and ending near the false vacuum: 

4>(0) f» 0tv , 0(o- ma x) ~ <j>tv ■ (2.13) 
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Other bounces exist [12J, [l3j and can be treated with completely analogous methods. Bounces can be labeled by an 
integer n characterizing how many times they cross the barrier. We will refer to the Coleman-De Luccia bounce 
as a "single bounce" solution, and the n > 2 bounces will be termed "oscillating bounces" 0|. I n the flat space 
limit, only single bounce solutions have finite action. The explicit bounce solutions can be found numerically by a 
straightforward shooting technique, as follows. 
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FIG. 2: Single bounce solution to 12. 31 - 12. 61 . for parameter values: j3 = 45, b = 0.25, e = 0.23. The shooting procedure 
determines the initial scalar field value (pi(0) = tpty — 0.000145 223 243 267 576 9, and u max is given by: HtopCmax = 3.7475. 
Here sim = sin (Hia), and Hi = 7r/cr^ ax . 




FIG. 3: Triple bounce solution to 12. 3H - 12.61 . for parameter values: f3 = 45, b — 0.25, e = 0.23. The shooting procedure 
determines the initial scalar field value ip 3 {0) = </Hv - 0.008 860 088 903 713 227, Htopffmax = 3.1243. Here sin 3 = -±- sin (H 3 a), 
and Umax is given by: H 3 = TY/a^^. 



The shooting parameter is the initial value </>(0) of the scalar field. This value is chosen near the true vacuum 
value and adjusted until the coupled initial value problem l|2.3|) - l|2.6|) has a solution satisfying both a(cr max ) = 




FIG. 4: Quintuple bounce solution to (I2.3H - 12.6L for parameter values: (3 = 45, b — 0.25, e = 0.23. The shooting procedure 
determines the initial scalar field value (ps(0) = <ptv — 0.337102 748195 2641, HtopcSuc = 3.1264. Here sins = ^ sin (H^a), 

and iTmax is given by: H 5 = it/a max . 



and </>(<r max ) = 0, for some cr max . The value of cr max is determined by this shooting procedure, and so depends on the 
bounce. Since the metric field behaves as a(er) ~ a for small cr, we cannot directly start integrating (|2.3J) at a = 0. 
Instead, we Taylor expand both a(a) and <f>{cr) about 0, and use these Taylor expansions to begin the integration at 
a point very close to 0. We then do a shooting scan of (f)(0), adjusting it digit by digit, in a rational form to preserve 
precision. This computation is simple to implement in Mathematica. In a few minutes one can determine </>(0) to 
32 decimal places. We found that the shooting went faster with the following simple rescaling of the differential 

equations, as in ^5|: we rescale ip = 4>/v, and a and a(a) as: s — J (3H^ op a, and p(s) = *J (3H^ op a(a), so that the 

parameters v and H top scale out of the classical equations of motion. 

Representative examples of bounce solutions with (j>(0) ~ tfttv and 0(er max ) ~ (f>i v are shown in Figures [21 El and 01 
These plots are for (3 = 45. For larger (3 more oscillations are possible in the bounce solutions. Note that the metric 
field a(a) deviates significantly from the sinusoidal form of l|2.12|l for the single and triple bounce, but less so for the 
quintuple bounce. This, together with the fact that the scalar field (f> is closer to the Hawking-Moss constant value of 
'/'top = 0, is another reflection of the fact that the highly oscillating bounce solutions tend towards the Hawking-Moss 
solution [Tflj . 



III. FLUCTUATION OPERATOR 



Having reviewed how to find the classical bounce solutions, we now turn to the problem of fluctuations about these 
solutions. With the inclusion of gravity, the fluctuation problem becomes more subtle, and the fluctuation operator 
about such cosmological instantons has been widely studied [2(| H^, |2^, |3(| ■ The variation of the action under 
fluctuations of the scalar field and the metric field requires a nontrivial constraint analysis, with different possible 
gauge fixing procedures. For the purpose of this paper, we choose the gauge fixing scheme described in Section 4.3 
of [29|. and Section IV of [23]. In this scheme, the only physical degree of freedom is the fluctuation of the scalar 
field, and the second variation of the action can be expressed as (compare with Eqn. (12) in j2]j and Eqn. (18) in the 
second reference of H3): 

S^M= 7 r 2 j d CT ^|-A^^M^ + a 3 ((T)f/[a(a))0(a) ^|^ . ( 3. 1} 



Here 



and 



r\( T ) — i , *°W(*) f „ 2 v 

s(CT) = 1 + 2(-A 3 -3iq ' (3 - 2) 

U[a,4>] = -V W + — + K {— + Q2( -A 3 -3if) ' (3 - 3) 
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Here we have temporarily re-instated the K dependence, although in our numerical studies we return to K = 1, and 
A 3 is the Laplacian on S 3 . 

To pass from this secondary action to the Jacobi equation j^, a Sturm-Liouville differential equation whose 
eigenvalues will define the determinant of the fluctuation operator, we need to specify the weight function. The 
weight function is determined by defining 

= J d 4 x^g(S(t)) 2 = 2tt 2 J daa 3 (a) [8<j)(a)] 2 (3.4) 

Then in terms of the perturbation function $ = 5<j), the fluctuation equation (the Jacobi equation is 

1 d /a 3 d$ 
a 3 da \ Q da 



73- (tt— )+U[a,ct>}$ = \Z , (3.5) 



which is defined on the interval a G [0, cr max ], with Dirichlet boundary conditions, and where A denotes the eigenvalue. 
The "fluctuation potential" U[a,cj>] is the function in 1)3.3(1 . The S 3 Laplacian A3 appearing in 1(3.2(1 and 1(3.3(1 can be 
replaced by its eigenvalue: —A3 — > 1(1 + 2), so we obtain a fluctuation equation as an ordinary differential equation 
for each integer value of I (Z 7^ 1). 

In the flat space limit, k — > 0, Q — ► 1, a — ► a, and cr m ax — > 00; in which case we recover the familiar flat space 
fluctuation equation [Tflj 



r 3 dr \ dr 1 



U[(/>(r)]<f> = A$ 



U[<t>) = V"{<t>)+ 1 -^- , (3.6) 

where a becomes identified with the Euclidean length r, which ranges from to 00. Much is known about solutions 
to this flat space fluctuation equation 1(3.6(1 . Our goal now is to study some properties of the more general fluctuation 
equation in ((3.5|) . 

For completeness, we note here that for the purposes of discussing the existence of negative modes it is possible to 
make other choices of the weight function, which yield superficially different-looking Jacobi operators. In Appendix A 
we give the explicit transformation between our choice 1(3.4(1 of weight function and those made in [2^, |27], |28l |2^, 03 • 



IV. LARGE I BEHAVIOR OF FLUCTUATION DETERMINANTS 

Both <p(a) and a(a) are functions just of the proper time a, so the fluctuation problem separates into partial waves, 
which can be labeled by an integer Z, just as in the flat space case. Then formally we can write the log determinant 
of the fluctuation operator A as 



In 



Det [A] 

Dct [Afree 



]T(Z + l) 2 ln 



Det [A«] 



Dct 



A 



(0 

free 



(formal !) , 



(4.1) 



where A^) is the differential operator in 1(3.5(1 . for each Z, and the eigenvalue of — A3 is 1(1 + 2), with degeneracy (1 + 1) 2 . 
This formal expression 1(4.1(1 must be interpreted with caution, because as in the flat space case, for low Z values there 
may be negative and zero modes. Nevertheless, for generic Z, since each A^ is a one-dimensional differential operator, 
the determinant ratio ^Det [A'''] /Dct A^ ^ is finite, and can be computed efficiently using the Gcrfand-Yaglom 
technique |33L 13^ . I35L I36L l37j . In this approach one simply numerically integrates both Jacobi equations, A^ $W = 



- 0, for zero eigenvalue and with suitable common initial value boundary conditions, <&">(a) 
(7 = 0. Then the ratio of the determinants is the ratio of these two functions evaluated at a mayi : 



and A^ rc 



Dct [A 



(0 



Dct 



A 



free 



$L(<Xmax) 



a 1 at 



(4.2) 



This technique provides a simple computational method for evaluating the finite determinant for each Z, without ever 
having to compute any eigenvalues. However, of course, even though each term on the RHS of 1(4.1(1 is finite, the 
sum over Z diverges |35j. This is clearly because we have not regularized and renormalized the determinant. This 
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divergence is not a feature of the gravitational coupling — exactly the same thing happens in flat space |23|, [24| , 
where one can indeed extract a finite renormalized determinant by subtracting certain known contributions from 



In (Det [A^] /Det 



A 



(0 



free 



for each I, rendering the sum finite. The precise form of the subtractions can be found 

in various ways, using Feynman diagram techniques p3[ . zeta function regularization [36|, or radial WKB 123. The 
finite part of these subtractions is related to the specific renormalization prescription |22l I23I I24I l38l l39l | . In |24j , in 
flat space, it was checked explicitly that the result of this procedure connects smoothly to the analytic thin- wall limit 
results for the renormalized fluctuation determinant. 



A key element of this approach is knowledge of the large I behavior of In I Det [A^] /Det 



A 



(/) 

free 



which must be 



subtracted to make the I sum finite (renormalization involves a further step). In flat space the radial WKB analysis 
leads to the following expression for the large I behavior |24j : 



In 



Det [A 



Det A 




i / °° dr r W | f~ dr r 3 W{W + 2V" [fc] ) 



{i + iy 



o 



(i + iy 



(4.3) 



Here W = V"[(/>] — V'[^>f v ]. Subtracting these terms makes the sum over I in (|4.1|l finite, and is one part of the 
analysis leading to a finite renormalized fluctuation determinant. We now turn our attention to the large I behavior 
of these determinants in the gravitational case. 

It is immediately clear that with gravitational coupling such a computation cannot be done for the renormalized 
fluctuation determinant, as we do not know how to renormalize gravity. Nevertheless, we can study this question in 
the de Sitter limit, where the gravitational background is fixed to be of the de Sitter form in (|2.12() . This limit is 
physically appropriate when the variation of the potential on the scale of the barrier is much less than Vt op |4ll | . For 
the moderate values of the cubic coupling b in (|2.7|) considered here, this amounts to the condition /3e 2 <C 1, which 
means that the potential is large and positive, with: 



fv 



^top * H t , 



H 



(4.4) 



Thus, fixing the metric to have the de Sitter form 

a(cr) 



H 



sin {Ha) 



the classical equations of motion reduce to a single equation for 4>{a): 



(4.5) 



4> + 3~j>-V'(<t>) =0 
a 



(4.6) 



Even though a{a) is determined, one still finds various different types of oscillating bounce solutions for (f>{a), as 
described in the previous section. These have been extensively studied recently in [T^ | with a different choice of 
parameter /? = 70.03. 

To compute the determinant ratio in Q4.2[l . consider first the false vacuum case, which is the appropriate "free" 
reference operator: A^ o = AfJ. Since <j> = 0f v is constant, and V'{<pi v ) vanishes, the fluctuation potential 1)3. 3|) 
simplifies dramatically, and the Jacobi equation (|3.5|l for zero eigenvalue becomes 



V"{cp iv ) 



l{l 



4> 



(/) 

fv 



(4.7) 



j_ 

fffv 



sin{Hf v a). The zero mode solution with the correct initial value behavior, {a) ~ a 1 , is an 



with a{a) 

associated Lcgendrc function (essentially derivatives of a conical function) 



A^fv 



P 



i+i 



sin {H fv a) 



[cos {H iv a) 



(4.8) 



where Nf v is an unimportant normalization constant. This function is positive definite and diverges as Hf v a — > ir. 

Now consider the same computation but for a bounce solution. Immediately we find a significant difference between 
the flat and gravitational cases. In flat space both free and bounce solutions are defined on the same interval r £ [0, oo). 
But with gravity, a nontrivial bounce solution is defined on the interval [0, cr^™ nce ], where the interval is determined 
by the second zero of the metric function a{a). So, in general, for solutions of the full bounce equations (I2.3|) - (|2.5|l . 
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1 

FIG. 5: In [DetA (0 /DetA^ e ] for the single bounce, for parameter values: = 45, b = 0.25, e = 0.046, ip(0) = <p tv - 

0.000 551736 739 278 415 5, <p(ir) = tptv + 0. 000 1 929 113 943 772 119 407 698 5. Diamond points denote the numer ical r esults 
using the Gel'fand-Yaglom technique as in 14. 1011 - 14. 1 IB : the (blue) dash line denotes the leading large I behavior in 14. 1211 with 
a defined in 14.141 and 7 = 0; the (red) solid line denotes the leading and subleading large / behavior in 14. 121 with ot in 14. 141 
and 7 defined in 14.151 . 



the false vacuum solution and a nontrivial bounce are defined on different intervals. Fortunately, this problem goes 
away precisely in the de Sitter limit being considered here, where we can take the metric field to be of the form in 
14.51) with H = Hf v , so that both solutions live on the same interval. 

As in the flat space case [23|, |24| , given that the free solution (|4.8|l is known analytically, it is better to consider the 
ratio of the functions appearing in (|4.2I) : 



= $('>(<r) 



(4.9) 



because this ratio remains finite as a 
value boundary conditions: 



Cmax- This ratio satisfies the following differential equation with simple initial 



T 



fv 



fv 



T 



QU 
T(0) = 1 



Q £tv 
Q $fv 



U {v T 
T(0) = 



(4.10) 



This also shows why the normalization of the false vacuum solution $f v in l|4.8|l is not important. It is conventional 
to compute the logarithm of the determinant ratio, in which case the Gel'fand-Yaglom result (|4.2() can be written 
simply as 



In 



Det [A«] 



Det A 



(/) 



In 



r«(a max ) 



(4.11) 



We have computed these determinants as a function of I, by numerically integrating the initial value problem (|4.10f) 
for various bounce solutions, using the fluctuation equation in (|3.5f) . with metric field given by the de Sitter form 
14. 5|) . and the scalar field given by solutions to the bounce equation l|4.6|l . The scalar field bounce solutions in this 
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FIG. 6: In |DetA (0 /DetA^ e I for the triple bounce, for parameter values: (3 = 45, b = 0.25, e = 0.046, <p(0) = ip tv - 

0.007 517 937 804 678 529, <p(ir) = ip iv + 0.00 1906 935 294 979 618. Diamond points denote the numeri cal re sults using the 
Gel'fand-Yaglom technique as in 14.101 - 14.111 : the (blue) dash line denotes the leading large I behavior in 14.121 with a defined 
in 14.141 and 7 = 0; the (red) solid line denotes the leading and subleading large / behavior in 14.121 with a in 14.141 and 7 
defined in 14.151 . 



de Sitter case with same parameter (3 = 45 are very close to those shown in Figures [3 [3] and 0] so we do not bother 
plotting them again. The results for the logarithm of the determinant ratio are shown as diamond points in Figures 
[5] through [SJ The first empirical observation is that the large I behavior is very similar to that l|4.3[l found in the flat 
space case [24|: 



In 



Det(A«) 
Det(A[l) 



l + l (l + lf 



00 



(4-12) 



To extract the leading coefficient a, we consider the leading and subleading large / behavior of the Jacobi equation 
(13311 : 



_$(0 -3-6W 



V"{4>) + 2k<P 2 



1(1 



$(0 = A G)$« 



(I » 1) 



(4.13) 



Adapting the WKB analysis of [2J, |38j leads to the following result for the leading large I dependence of the log of 
the determinant ratio: 



daa{a) V" (</>) + 2k4> 



12 



V" 



9fv - 



(4.14) 



Notice the close similarity to the leading term in the flat space large I behavior in l|4.3|) . This curved-space leading 
behavior a/(l+ 1), with a given by l|4.14|l . is shown in Figures through |S] as dashed (blue) curves, and we see that 
the agreement at large I is very good, ft is much harder to find the next-to-leading behavior because the subleading 
I dependence of the Jacobi equation is very complicated. Nevertheless, by analogy with the flat space case 1)4. 3J) we 
propose the estimate 



da a 3 (a) V"{<j>) + 2 K (j) 2 - V"((f) iv ) V"(4>) + 2n<f) 2 + V"{<j> fv ) 



(4.15) 
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FIG. 7: In petA w /DetA^ c J for the quintuple bounce, for parameter values: j3 = 45, b = 0.25, e = 0.046, ip(0) = tp tv - 

0.303 787 546 677 413 1, v p ( 7r ) — ftv + 0.187 525 348 746 153 5. Diamond points denote the numerical results using the Gel'fand- 
Yaglom technique as in 14.10H - 14.llll : the (blue) dash line denotes the leading large I behavior in 14.121 with a defined in 14.14H 
and 7 = 0; the (red) solid line denotes the leading and subleading large / behavior in 14.1211 with a in 14.141 and 7 defined in 

Including this subleading behavior in (|4.12|) produces the solid (red) curves in Figures [S] through [SJ and we see that 
the agreement with the exact numerical results is noticeably improved and is excellent at large I, and is characteristic 
of an asymptotic large / expansion in its behavior at small I. Thus, the estimate in l|4.15|l is very close to the exact 
answer. We found similar behavior for other bounces. 

V. NEGATIVE MODES 

In this section we turn to another important property of the fluctuation operator, namely the existence of negative 
modes. Here, to be more general we can return to the general bounce solutions, not needing to work in the de 
Sitter limit any more (although we find the results to be the same in either case). In the flat space false vacuum 
decay problem, it has been shown that the fluctuation problem <|3.6|) has one and only one negative mode, and that 
this occurs in the I = sector This single negative mode plays an important physical role in the semiclassical 
quantization, accounting for the decaying nature of the process [l&lisllisj . In the gravitational case, there has been 
considerable effort analyzing the appearance of negative modes in the / = sector |2(| I27L l28l l29t l30t l31| . For the 
scalar field fluctuations characterized by the secondary action l|3.1|l . the oscillating n-bounce solution has n negative 
modes in the I = sector j^J. Here we show that the higher n oscillating bounce solutions also can have negative 
modes in higher I sectors, while the single bounce solution, the Coleman-Dc Luccia solution, has precisely one negative 
mode, which is for I = 0. 

A direct numerical method for counting negative modes is based on an important theorem in the calculus of 
variations due to Morse |32l | . It states that the number of negative modes of the fluctuation operator A is given by 
the number of zeros of the solution of the initial value Jacobi equation A$ = 0. This is consistent with the related 
GePfand-Yaglom result (|4.2|l for the computing the determinant as the value of 4 > (tT ma x), since an odd number of 
zeros leads to a negative determinant. This Morse analysis has been applied to the counting of negative modes in the 
flat space false vacuum decay problem in [43j . 

So to count the number of negative modes for a given bounce solution, we numerically integrate the fluctuation 
Jacobi equation (|3.5[) . with initial value boundary condition <f> ( - i ' ~ a 1 , and count how many times this function 
changes sign on the interval [0, cr max ]. In this case we can do this computation using the full bounce solutions as 
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FIG. 8: In [DetA (i) /DetA^J for the Hawking-Moss bounce, for parameter values: j3 = 45, b = 0.25, e = 0.046. Diamond 

points denote the numerical result using the Gel'fand-Yaglom technique as in (I4.10H - 14.11II : the (blue) dash line denotes the 
leading large I behavior in 14.121 with a defined in 14.141 and 7 = 0; the (red) solid line denotes the leading and subleading 
large I behavior in 14.121 with a in 14.141 and 7 defined in 14.151 . 



TABLE I: The number of negative modes without including the degeneracy factor (I + l) 2 , for the parameter values: (3 — 45, 
e = 0.23. 
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obtained in Section^ not just in the de Sitter limit where the metric function a(a) is chosen to take the sinusoidal 
form (|4.5J) . In this way, we confirm the results of [27Ll3ll] that the n-bounce solution has n negative modes in the / = 
sector. More surprisingly, we find that for I > 2 there are some negative modes for the higher n oscillating bounce 
solutions. The precise pattern depends on the parameters in the potential, but a representative counting is shown in 
Table [I] As the oscillation number of the bounce increases there are more negative modes, and they extend to higher 
values of I. In studying many single-bounce solutions we have always found only one negative mode, and always in the 
/ = sector. We also found the same negative mode counting pattern using Lavrelashvili et aVs fluctuation operator 
in (Q)l . 

To put this result of extra negative modes at hi ghe r I in some perspective, consider the Hawking-Moss solution, 
which is the large n limit of the n-bounce solution Here, we can write the exact solution to the zero eigenvalue 
Jacobi equation, A<1> = 0, with initial behavior $ ~ a 1 , as an associated Legendre function, analogous to the false 



[cos {H top a)\ . (5.1) 



vacuum solution in 



Here /3 is the parameter defined in (|2.10(l . and N-rm is an unimportant normalization constant. The counting of the 
zeros of this function can be done precisely, and one finds that the number of zeros depends critically on the value of 
/3. For N(N + 3) < (3 < (N + l)(iV + 4), there are N +1 zeros. By the Morse theorem the number of nodes of this zero 
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mode solution is equal to the number of negative modes in the perturbation. This counting is also shown as the last 
row in Tabled and we have also confirmed that the numerical integration and the exact result give the same counting. 
Since the oscillating bounce solutions tend to this Hawking-Moss solution, this goes some way towards explaining the 
origin of these new negative modes for the oscillating bounce solutions at higher I. Physically, this is extra evidence 
that these higher n oscillating bounce solutions are not directly related to quantum tunneling, as suggested already 
in [13. 

VI. CONCLUSION 

In this paper we have analyzed several issues concerning the quantum fluctuations about classical bounce solutions in 
the theory of a self-interacting scalar field interacting with gravity. In flat space the semiclassical fluctuation analysis 
can be done completely, both analytically in the thin-wall limit and numerically for more general potentials. In the 
gravitational case, the fluctuation problem still separates into a set of one-dimensional fluctuation problems labeled by 
an integer I. We found the leading large I behavior l|4.12|l - (|4.14|l . and an estimate l|4.15|) for the subleading behavior, 
of the logarithm of the determinant of the fluctuation operator. The agreement with the numerical computations 
is impressive. We also analyzed the existence of negative modes using Morse's theorem, confirming that the single- 
bounce Coleman-De Luccia solution has a single negative mode, which lies in the / = sector, and that the oscillating 
n-bounce solution has n negative modes in the I = sector. We also found new negative modes for the oscilla ting 
n-bounce solutions for higher n with I > 2. This adds further weight to the physical interpretation suggested in [pj 
that these bounces are not directly related to quantum tunneling, but rather are related to the thermal character of 
quantum field theory in de Sitter space, and interpolate to the Hawkin g-M oss solution for large n. 

Many problems remain. The standard scalar fluctuation analysis [26|, |271 |2& u3 m the gravitational case precludes 
consideration of the I = 1 sector, and so we cannot yet say anything about the collective coordinate contribution to 
the renormalizcd fluctuation determinant. In the flat space case it was recently shown how this I = 1 contribution, 
combining the determinant with the zero modes removed and the collective coordinate contribution, could be expressed 
simply in terms of the asymptotic properties of the classical bounce solution j24|. Whether something like this can 
be found for the gravitational case depends on a different analysis of the I = 1 fluctuation problem. Perhaps the 
most challenging problem is that the renormalization of quantum gravity is not understood. In the flat space case, 
without gravity, the subtractions made from the regularized determinant for each I include a finite piece depending 
on the regularization scale. These can be associated with renormalization, permitting the computation of a finite and 
renormalized fluctuation determinant |23l l24l | . An important preliminary step for the gravitational case would be to 
develop fully this approach in the limit where the gravitational background is fixed to be de Sitter, in which case the 
large I behavior of the log determinants is given by l|4.12|l . and where the pcrturbativc renormalization of the scalar 
field in a fixed curved background is known [icj. Hopefully this can shed further light on the important question of 
the nature of the semiclassical path integral approximation in the presence of de Sitter gravity [U |4Ji US • 

Acknowledgments: Wc thank the US DOE for support through the grant DE-FG02-92ER40716. 

Appendix A . RELATED FORMS OF THE FLUCTUATION OPERATOR 

In discussing the existence of negative modes it is possible to make other choices than (|3.4ll for the weight function, 
yielding superficially different-looking Jacobi equations [2^, [2?], H^, 12^. l3l| . But for the purposes of computing the 
determinant, where the magnitude of the eigenvalues is also relevant, the choice in (|3.5|) is the most direct. For 
completeness, the choice of Lavrelashvili et al is to use the perturbation function = \J a 3 /Q5<j), and weight 
function a 3 /Q, leading to the Jacobi equation psl l2flL Isij 



d 2 $ 



L 



da 2 



U L [a,<f>}$ L = \ L $L , (A.l) 



where 



1 // , ^ Sua 2 , ., % ,9 &naa(b ,, . k 

UlM , -yea - 2Q2( _ A3 _ m { vw + Q2( _ A3 * m n*) - - 6 ^ 2 +vw 



1 10 12 \ 1 



(-A 3 -3K)(Q + 2)- - 2A3 - 8K 



(A.2) 
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which agrees with Eqn. (19) in the second reference in [3l| when K = 1 and A3 = 0. Furthermore, this potential 
satisfies 

U L [a,<t>] = QU{a,ct>] + ^-^^ . (A.3) 

Using this fluctuation operator, we found the same pattern of negative modes shown in Tabled 

Another choice, yielding an elegant result for the existence of negative modes, was made by Turok et al, who chose 
the perturbation function $t = Scj}/^ 1 , and weight function l/</> 2 , leading to the Jacobi equation [2rl \2l\ 

~(^-^j+U T [a,^ T = XT^T , (A.4) 

Here 

U T [a,<t>] =a<p 2 {-A 3 -3K) , (A.5) 

and it is related to our U[a, (f>\ by 

C/ T K0] = a 3 2 C/[a,0]-^ f^J . (A.6) 
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